Investigation of exciton transport in crystalline thin-films of the organic 
semiconductor di-indeno-perylene using photoluminescence analyses 
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The exciton transport in the prototypical organic semiconductor Di-Indeno-Perylene (DIP) has 
been investigated by means of photoluminescence (PL) quenching and interpreted by an advanced 
exciton diffusion model including interference effects, quencher penetration as well as non-perfect 
exciton quenching. X-ray diffraction revealed a correlation between the exciton diffusion length of 
about 90 nm and the structural coherence length of the DIP layers. Temperature dependent studies 
in a range of 5 - 300 K indicated an incoherent exciton transport above 80 K at activation energies of 
10 - 20 meV related to the thickness dependent gradient of exciton density. Below 80 K a coherent 
exciton transport can be observed by the reduced phonon-interaction at cryogenic temperatures. 
This manuscript is a pre-final version. 



By their progressively improving power conversion ef- 
ficiency, organic thin film photovoltaic cells pave their 
way towards first applications^. In organic bilayer pho- 
tovoltaic cells, the competition between exciton diffusion 
and absorption length is decisive for the applicable thick- 
ness of the photo-active layers and thereby the overall 
device efficiency^. As it became obvious from observa- 
tions of band-like, coherent exciton motion in organic sin- 
gle crystals^ that morphology significantly influences the 
diffusion characteristics of optical excitations, molecular 
model systems supporting coherent exciton transport by 
their crystalline film structure are highly desired 87 . To 
aim for this challenge we have performed temperature 
dependent exciton transport studies in thin films of the 
prototypical organic semiconductor Di-Indeno-Perylene 
(DIP). 

Bilayer thin films cells of DIP in combination with the 
electron acceptor Ceo have shown superior performance 
including fill factors above 70 % and power conversion 
efficiencies of up to 4 %2L. 

Therefore, in this paper we present an approach, 
including morphological aspects as boundary effects, 
quencher penetration and optical interference effects, to 
determine the excitonic transport properties by thick- 
ness dependent Photoluminescence (PL) studies. This 
method also offers insight in the dominate excitonic 
transport mechanisms by temperature dependent PL- 



studies. For this purpose, DIP samples were prepared by 
organic molecular beam epitaxy on cleaned glass slides. 
The resulting film structure of the samples was ana- 
lyzed in detail by X-ray diffraction at various geometries. 
Bragg-Brentano scans confirmed the previously reported 
(001) lattice spacing of about 1.66 nm corresponding to 
an almost upright orientation of the moleculesA and the 
long-range coherence of the lattice planes is indicated by 
pronounced Laue oscillations at the (100) Bragg peak 
(left inset of fig. [TJd). Furthermore, DIP film mosaic- 
ity is charaterized by an averaged crystallite tilting of 
only 0.04 deg along the surface normal, i.e. along the 
exciton transport direction (right inset of fig. [TJd). To 
define a measure for the film quality we employed the 
ratio of specular versus diffuse intensity (dark vs. light 
blue shaded areas in the right inset of fig. Q~|d), the lat- 
ter being influenced by structural inhomogeneities within 
the organic film 13 . As can be seen, this ratio peaks at 
a thickness of about 70 nm in agreement with the crys- 
tallite height, determined by Laue-oscillations. Based on 
this data a model of the underlying film structure has 
been deduced and is described by a complementary error 
function representing the integrated Gaussian distribu- 
tion of DIP crystallite heights (fig. [TJ;). 

For the Photoluminescence (PL) studies DIP thin films 
were prepared with Copper Phthalocyanine (CuPc) ther- 
mally evaporated on one half of the DIP stepped wedge 
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FIG. 1. (a) Structural model of DIP. (b) X-ray structural 
analyses and Laue-oscillations (upper left inset) indicating 
the high crystallinity of the DIP films. The ratio between 
specular intensity (dark blue area) and the diffuse scattering 
intenity (light blue area) has been considered as a measure 
for the structural quality of the DIP layers (upper right in- 
set), (c) Model of the DIP film texture together with possible 
exciton recombination processes. Integration of the Gaussian 
distribution of the DIP crystallite heights (black curve) yields 
the complementary error function (red curve). 



(inset of fig. [5]), acting as exciton quencher due to its 
HOMO and LUMO positions of 5.1 eV and 3.2 eV with 
respect to those of DIP, 5.4 eV and 2.9 eV 7 ' 21 . The PL 
intensity of the DIP films, illuminated through the sub- 
strate by a NdYag-Laser at 532nm, was measured by a 
sensitive CCD detector after passing a 590 nm long pass 
filter to block background intensity of the excitation. For 
a given DIP layer thickness, the comparison of the rela- 
tive PL intensity ratio Q=PL® jPL n ® , recorded on the 
bare {PL 71 ®) and the CuPc-covered (PL®) area reveals 
contributions from non-radiative decay channels and dis- 
sociation of excitonsi 7 -. In films thinner than Lo most 



of the excitons will reach the quenching interface within 
their radiative lifetime, thus strongly reducing the PL 
signal. On the other hand for films thicker than Lo the 
PL signal is expected to dependent only on the absorp- 
tion and therefore should be identical for the free and the 
quencher- covered areas. As a result, Lo in the molecu- 
lar layer can be determined by thickness dependent PL 
measurements. It is important to stress at this point that 
according to the layer thicknesses and Lo under studyi& 
as well as by the small spectral overlap for DIP, contri- 
butions to quenching by Fluorescence Resonance Energy 
Transfer (FRET) can be neglected^. 

However, the situation becomes more complex for real 
samples, especially if the absorption layer thickness ap- 
proaches the excitation wavelength within the organic 
material. For thicknesses equal or larger than nA/4 in- 
terference effects have to be taken into account 2 ^, requir- 
ing detailed information on the indicatrix of the organic 
material. Therefore, to quantify exciton quenching one 
has to calculate at first the intensity profile of incident 
light across the sample yielding the initial exciton density. 
Accounting for light reflection at the organic interfaces 
yields the depth dependent intensity profil o 19 ' 20 

= e ~ az + p 2 e~ a(2d ~ z) 

+ 2pe- ad cos(-2k(d - z) - 5) (1) 

Here z describes the position in the film (z = refers 
to the glass/organic- interface), a is the absorption coef- 
ficient, p is the reflection coefficient at the organic in- 
terface, d is the sample thickness, k is the wave vector 
and 5 is the phase shift occurring upon reflection at the 
backside of the organic layer. With n(z,t) being the local 
exciton density at time t, the continuity equation reads 



^,t) ^ D d 2 n{z,t) _n^ +GM (2) 
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For given exciton lifetime t, the exciton generation 
term G(z,t) under constant illumination, i.e. at steady 
state conditions (n(z, t)/t = 0), is expected to be propor- 
tional to the light intensity profile. Appropriate bound- 
ary conditions can be achieved by neglecting quenching 

and as- 
= 0, or no 



at the glass/organic- interface, 
suming either complete quenching, n(z,t) 
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dz 
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quenching, 
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0, at the organic/quencher- or 
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organic/air- interface, respectively. Real samples will al- 
ways operate between these two extremes, i.e. partial 
quenching might occur at the organic/air- interface e.g. 
due to surface states as well as non-ideal quenching at 
the organic/quencher- interface. For perfect quenching 
the absolute value of the slope and accordingly the ex- 
citon concentration gradient will be at maximum S max 
at the quencher interface. This provides a measure for 
the quenching quality V = S/S ma x, where S is the 
slope at the organic/quencher- interface assuming incom- 
plete quenching. Obviously, V assumes values between 
and 1. As for more complex intensity profiles eq.@ 
can no longer be solved analytically, we employed an 
ansatz based on superposition of Dirac delta functions 
G{z) = 5{z — zq)/t as generation term. The low ex- 
citon densities in our PL studies render the diffusion to 
be concentration independent and therefore n$(z, zO) can 
be used in the generation profile. A convolution of the 
light intensity profile, including interference effects and 
rig(z, zO) results in 

d 

n(z) = J $ • I(zq) ■ n s (z, z )z (3) 



The factor $ describes the quantum efficiency of exci- 
ton generation upon photon absorption. Assuming non- 
radiative decay processes in the absorption layer to be of 
minor importance, the detected PL intensity corresponds 
to the number of excitons not quenched within their ra- 
diative lifetime 



Here we used the common definition of the exciton dif- 
fusion length Ld = V Dr. For V = each exciton con- 
tributes to photoluminescence thereby yielding the inte- 
grated intensity profile. Using eq.(jl]) the PL signal can 
be calculated for any given I{zq). The thickness depen- 
dent Q measured at room temperature (RT) is displayed 
in fig. [5J Three distinct features can be clearly identified: 
First, Q does not approach zero even for doip — > 0, sug- 
gesting imperfect quenching at the DIP/CuPc interface. 
Second, the oscillatory thickness behavior of Q refers to 
the presence of interference effects. Third, even for film 



thicknesses of 300 nm Q is still below 0.6, indicating the 
existence of a significant, almost thickness independent 
quenching contribution. 
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FIG. 2. Thickness dependent quenching ratio 

Q=PL Q /PL nQ . The fit (blue line) corresponds to the 
model discussed in text accounting for interference effects 
and interdiffusion at the DIP/quencher interface. The inset 
displays the sample structure. 

To relate these results to the underlying film morphol- 
ogy, complementary studies on the temperature depen- 
dent PL have been carried out on DIP layers of vari- 
ous thicknesses. Fig. [3^ shows the PL intensity as a 
function of the inverse thermal energy recorded on thin 
(30 nm) and thick (290 nm) DIP layers with an with- 
out quenching layer. Below 0.145 (keV) -1 , i.e. above 
80 K, the fluorescence of the free (PL n Q) and of the 
CuPc covered (PL®) DIP volume is thermally activated, 
with PL n< ^ showing a distinct increase with thickness 
from 12 meV to about 21 meV. This tendency can be 
related to different exciton density gradients as a func- 
tion of layer thickness. A steeper concentration gradi- 
ent, occurring for thinner films, reduces the effective en- 
ergy required for transport^. In addition, this is con- 
firmed by the quencher covered 30 nm thick DIP film 
whose PL® activation energy is about four times smaller 
than that of the uncovered film. Finally, the value of 21 
meV for the thickest DIP film is in good agreement with 
the PL activation energy of covered DIP irrespective of 
layer thickness and, therefore, refers to exciton transport 
within the DIP crystallites. In analogy to charge car- 
rier transport studies on long-range ordered molecular 
stacks 5 wc attribute the temperature dependent fhiores- 
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cence behavior to exciton-phonon-interaction constitut- 
ing a non-radiative decay channel which is suppressed 
towards lower temperatures, i.e. upon freezing-out non- 
local phonon modest. This assumption is corroborated 
by the observed activation energy range which coincides 
with that of non-local phonon modes in polyaromatic sin- 
gle crystal a 5 i 16 . Furthermore, data for all thicknesses in 
fig. [3] reveals the PL intensity to become temperature 
independent below 80 K. Caused by an almost coherent 
transport of the (free) excitonic species within the crys- 
talline DIP fraction^ without dissipative effects by static 
inhomogenities . 
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FIG. 3. (a) Comparison of the PL(T) behavior with and with- 
out quenching layer for the thickest and thinnest film of this 
study, (b) Temperature dependent Q for five representative 
DIP layer thicknesses normalized to the room temperature 
value. Thinner layer show the signature of a thermally ac- 
tivated Q above 80 K whereas below this temperature, the 
ratio remains constant. 

As it becomes obvious from fig. [3j}, the thickness de- 
pendent Q is characterized by a temperature dependent 
decrease for thin DIP films followed by an temperature 
independent regime below 80K. In contrast, thick films 



(d > lOOnm) show an almost constant Q over the en- 
tire temperature range. This points directly at the mor- 
phological origin of this phenomenon which can be in- 
terpreted by the previously discussed spatial extension 
of DIP crystallites along the exciton transport direction, 
i.e. normal to the (001) plane (fig[lj. For thick films 
(d > lOOrtm) non-radiative as well as radiative recom- 
bination processes take place during exciton transport 
within a single crystalline domain, due to a substan- 
tial energy barrier at the grain boundaries. However, 
for thin films (d < 100 wn), where most of the crystal- 
lites will stretch across the entire film, exciton transport 
and recombination are substantially determined by the 
respective quenching interface. Consequently, the vari- 
ation with temperature is mainly governed by the PL® 
contribution to Q. Furthermore this model is supported 
by the fact that the temperature independent Q below 
80 K is significantly smaller for thin than for thick DIP 
films as quenching is still present in case of the former 
whereas in thick films excitons are not being capable of 
reaching the quenching DIP/CuPc interface and there- 
fore relax to the ground state by the same mechanisms 
as for the uncovered DIP section. 

Finally, this model also explains the non-vanishing 
quenching contribution for very thick films (fig. [2]). 
In this case, the intrinsic crystallite size, which agrees 
very well with the optimized thickness deduced from the 
Ispec/Idiff ratio measured by X-ray (fig. Hh), initiates 
a pronounced 3D island growth (Volmer- Weber growth) 
together with an enhanced surface roughening as con- 
firmed by detailed x-ray analyses and atomic force mi- 
croscopy (AFM) 2 -. As indicated in fig. [Tb this topography 
promotes the penetration of the top-deposited molecular 
quenching material into the grain boundaries, leading to 
a significant contribution by lateral exciton diffusion and 
dissociation to the PL quenching of about 50 % even for 
thicknesses larger than 250 nm. In addition, this trend is 
supported by the better overlap of the electronic wave- 
function and the molecular transition dipoles within the 
(001) plane, presumably entailing the lateral components 
of the tensorial diffusion constant to be larger than the 
out-of-plane component^. A lower limit for the lateral 
Ld of several tens of nanometres has been reported which 
overall matches our observation 2 - 2 . However care has to 
be taken as those values have been determined in the 
presence of metallic nanostructures which might influ- 
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ence the exciton dynamics e.g. by plasmon mediated 
non-radiative recombination. 

Summarizing these results, we suggest a relation to ac- 
count for the morphological impact of real quencher in- 
terfaces, thereby resembling the aforementioned comple- 
mentary error function 

Q' = Q(R+(l-R)-erfc(^^\) (5) 

The average crystallite height and surface roughness 
are described by do and a whereas R qualitatively con- 
tains all effects which limit the relative quenching in case 
of thick layers. Modelling the thickness dependent Q by 
our extended exciton diffusion model, comprising effects 
by optical interference and quencher interface morphol- 
ogy, reveals the fit shown in fig. [5] and an Lrj of 90 nm 
for DIP. This value is in good agreement with previous 
reports using spectral photocurrent measurements inter- 
preted by a model of Feng and Ghosh 3, -. Comparing the 
estimated average crystallite height of do = 85 nm with 
L d renders the correlation between these properties to be 
evident. Furthermore, the adjusted reflection coefficient 
ratio between the free and the CuPc-covered DIP volume 
can be deduced by fitting the measured data points under 



consideration of eq.(JT|) and amounts to 1.4. Though influ- 
enced by the respective interface morphology, this result 
as well as the relative phase shift S almost matches ir/2 
which is theoretically expected for a standing-wave with 
open end (organic/air interface) versus that with fixed 
end (organic/quencher interface). Finally, the discussed 
transport model also renders the thickness dependent in- 
terface roughness which according to the rapid roughen- 
ing model observed for DIP scales as a—dP with /3=0.75^. 
In summary, we have investigated the exciton transport 
in long-range ordered molecular thin films by PL quench- 
ing, explicitly accounting for interference effects as well as 
structural properties at the respective quencher interface. 
Based on a diffusion model in combination with comple- 
mentary morphological analysis we were able to describe 
the thickness and temperature dependent PL quenching 
behavior of DIP. The resulting thermal transport behav- 
ior indicated coherent exciton transport in a tempera- 
ture regime below 80 K. The diffusion length of 90 nm 
is in striking agreement with the average DIP crystallite 
height and, with prospect on implementation in organic 
photovoltaics, indicates the decisive role of structural or- 
der on the exciton transport in molecular semiconduc- 
tors. 
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